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ABSTRACT 

We explore large-scale hydrodynamics of H II regions for various self-similar shock 
flows of a polytropic gas cloud under self-gravity and with quasi-spherical symmetry. 
We formulate cloud dynamics by invoking specific entropy conservation along stream- 
lines and obtain global self-similar "champagne flows" for a conventional polytropic 
gas with shocks as a subclass. Molecular cloud cores are ionized and heated to high 
temperatures after the onset of nuclear burning of a central protostar. We model subse- 
quent evolutionary processes in several ways and construct possible self-similar shock 
flow solutions. We may neglect the mass and gravity of the central protostar. The 
ionization and heating of the surrounding medium drive outflows in the inner cloud 
core and a shock travels outwards, leading to the so-called "champagne phase" with an 
expanding outer cloud envelope. Complementarily, we also consider the expansion of a 
central cavity around the centre. As the inner cloud expands plausibly due to powerful 
stellar winds, a cavity (i.e., 'void' or 'bubble') can be created around the centre, and 
when the cavity becomes sufficiently large, one may neglect the gravity of the central 
protostar. We thus present self-similar shock solutions for "champagne flows" with an 
expanding central void. We compare our solutions with isothermal solutions and find 
that the generalization to the polytropic regime brings about significant differences 
of the gas dynamics, especially for cases of n < 1, where rt is a key scaling index 
in the self-similar transformation. We also compare our global polytropic self-similar 
solutions with numerical simulations on the expansion of H II regions. We further ex- 
plore other possible dynamic evolutions of H II regions after the initiation of nuclear 
burning of the central protostar, for example asymptotic infiows or contractions far 
from the cloud centre and the ongoing infall around a central protostar. In particular, 
it is possible to use the downstream free-fall solution with shocks to describe the dy- 
namic evolution of H II regions shortly after the nascence of the central protostar. We 
also give an analysis on the invariant form of self-similar polytropic flows by ignoring 
self-gravity. 

Key words: H II regions — hydrodynamics — ISM: clouds — shock waves — star: 
formation — stars: winds, outflows 



1 INTRODUCTION 

The compact (scale sizes of ~ 0.1 — 1 pc) and ultracom- 
pact (UC; scale sizes < 0.15 pc) H II regions are asso - 
ciated with massive OB stars (e.g.. iHabing fc Israeli [l979h . 
The UC H II region stage ( < 10^ yr) stands for a sub- 
stantial part of the relatively short main-sequence lifetimes 
of OB stars. Several radio surveys (e.g.. Wood & Church- 
well 1989; Fish 1993; Kurtz, Churchwell & Wood 1994) ob- 
served expansions of luminous H II regions with shock signa- 
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tures. "Champagne flow" models fe.g.. |Tenorio-TaEle 


|l979|; 


iTenorio-Tade. York & Bodenheimeij 119791: YorkI 1986 


) sue- 



cessfuUy explain the expansion of H II regions by considering 
a protostar formed in a cloud core, photonionizing and heat- 
ing the cloud, as well as driving a shock that accelerates the 
ionized gas to expand rapidly. Observations tend to support 
"champagne flow" models, such as Lums den fc Hoare (1999) 
for U C H II regions G 29.96-0.02 and iBarriauh fc JoncasI 
(|2003) for Compact H II region Sh 2-158. Champagne flows 
in clouds of la r ger sc ales have also been identified, such 
as [Foster et al.l l|2007h for the dense Galactic H II region 
G84.9-I-0.5 and Maheswar et al. (2007) for the classical H II 
region S236 in the cluster of OB stars NGC 1893. 
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lTenorio-Tag"i3 (|l979l ) classified champagne flows into 
cases R and cases D. For the case R, the ionization front (IF) 
created by the emergence of the central massive protostar 
rapidly breaks out from the dense cloud and leaves the gas 
behind it fully ionized. For the case D, the ionization front 
is 'trapped' inside a cloud, and produces an expanding H II 
region within a cloud. In the formation phase of H II regions, 
whether the IF is R-type or D-type depends on the initial 
grain opacity, the io nizing flux, and the initial density and 
the s ize of a cloud l|Franco. Tenorio-Tagle &: Bodenheimerl 
I1990I ). In the expansion phase, H II regions with an ini- 
tial mass density profile p cx r"' and I > 3/2 are 'den- 
sity bounded', where p is the mass density and r is the 
radius (see, e.g., Osterbrock 1989; Franco et al. 1990). Q 
It is also possible that the D-type IF changes to a weak 
R-type IF when / > 3/2. In such cases, the fully ionized 
cloud begins to expand and an outgoing shock forms. This 
is referred to as the "champagne phase" (Bodenheimer et 
al. 1979). Notably if a shock front encounters a steep neg- 
ative density gradient, for example, the edge of a cloud, 
asymmetric "champagne fiows" may occur, as observed by 
Lumsden & Hoare (1999). According to the VLA survey 
by Wood & Churchw eU (1989), 16% o f H II regions bear 
cometary appearance. [Arthur fc Hoard (|2006l ) numerically 
simulated "cometary champagne fiows". If I < 3/2, the H 
II region is 'ionization bounded', e.g., t he ultr aviolet radi- 
ation is trapped within a finite radius l|Osterbro ck 1989). 
In such cases, the ionized region should expand as t^/C'^^'J^ 
driving a shock that would accelerate the ambient medium 
into a thin shell (Franco et al. 1990). In this paper, we focus 
on the champagne phase of a cloud assumed to be 'density 
bounded', implying I > 3/2. We also assume that the cloud 
is fully ionized shortly after the onset of nuclear burning of 
a ce ntral mass i ve pro tostar. 

IShu et al.l l|2002l ') (also Tsai & Hsu 1995) investigated 
isothermal "champagne fiows" under spherical symmetry in 
the self-similar framework. By neglecting the gravity of the 
central massive protostar and assuming that a cloud initially 
stays at rest and gets heated by the luminous massive pro- 
tostar to a uniform high temperature, one may obtain self- 
similar expansion solutions connected with isothermal out- 
fiows with shocks; this corresponds to a "champagne fiow" 
in a highly idealized setting. The initial isothermal mass 
density profile scales as p oc r~^. In addition, for molecu- 
lar clouds with other possible initial mass density profiles, 
IShu et all l|2002l ) also ignored the self-gravity of a cloud com- 
pletely and proposed another self-similar transformation, re- 
ferred to as the 'invariant form', and obtained solutions for 
cases with an initial mass density profile p ccr~' , where the 
power-law index / is not necessarily equal to 2. In general, it 
is not realistic to suppose molecular clouds to be isothermal 
in many astrophysical situations. One specific example of 
cloud temperature measurement of a UC H II region NGC 
6334F undergoing a "champagne fiow" reveals a conspicu- 
ous temperature grad ient from the centre to the edge (e.g., 
iDe Buizer et al.ll2002l ). 

Energy sources and plasma coolings in molecular clouds 



^ In Franco et al. (1990), the power-law index of the mass density 
profile is denoted by w instead of I as adopted here to avoid 
notational confusions. 



are not completely known. We then approximate the en- 
ergy equation by a general polytropic equation of state 
p = K{r, i) p^, where p is the thermal gas pressure, 7 is 
the polytropic index and the proportional coefficient K{r, t) 
(related to specific entropy) depends on radius r and time 
t in general. Setting global constant, the equation of 

state simply becomes a conventional polytropic one. By ad- 
justing 7, we may model various situations of H II regions in 
molecular clouds. For example, for 7=1 and a constant k, 
our solutions reduce to isothermal ones. Since "champagne 
fiows" in a polytropic molecular cloud has not be en studied, 
we wo uld generalize the isothermal analyses of IShu et al] 
(2002) and of Tsai & Hsu (1995) to a polytropic description 
of self-similar "champagne fiows" . We shall provide the ba- 
sic formulation with the most general polytropic equation of 
state (i.e., specific entropy conservation along streamlines; 
Wang & Lou 2008) and present global "champagne flow" 
solutions with shocks for a conventional polytropic gas. 

IShu et al.l (|2002l ) introduced the Bondi-Parker radius 
as a measure for the effective distance of the central gravity. 
The Bondi-Parker radius is defined by rep = GM», / {2a?) , 
where Af« is the mass of the central gravity source and a 
is the sound speed of the surrounding medium, which is a 
constant for an isothermal cloud. The mass originally resid- 
ing within a radius ro is dumped into the star during the 
star formation. After the star formation, as the surrounding 
cloud becomes much hotter with a higher sound speed, the 
Bondi-Parker radius becomes much less than ro. Therefore 
for the gas in r > ro, the gravity of the central massive star 
may be neglected. This reasoning naturally leads to the pos- 
sible existence of a cavity around the centre of a molecular 
cloud, which we refer to as 'void' or 'bubble'. At t = 0, the 
void boundary is at r = ro. For an expanding cloud, the 
central void also expands. Indeed, a stellar wind drives also 
a principal shock and is capable of sweeping the surround- 
ing ionized gas into an expanding shell. Wood & ChurchweU 
(1989) identified central cavities in the shell or cometary UC 
H II regions, which are thought to be supported by stellar 
wind and radiation pressures in their survey. Lumsden & 
Hoare (1999) suggested a "champagne fiow" surrounding a 
hot s tellar wind bubbl e to interpret observations of G 29.96- 
0.02. IComeronI (|l997l ) numerically simulated the dynamic 
evolution of wind-driven H II regions with strong density 
gradients and found that features of classical champagne 
model are not substantially changed, except that the com- 
pression of the swept-up matter would, rapidly and partic- 
ularly in densest cases, lead to the trapping of the IF and 
inhibit the champagne phase. Therefore, the dynamic evo- 
lution of void expansion represents an important physical 
aspect of "champagne fiows". Recently, Lou & Hu (2008, 
in preparation) explore self-similar solutions for voids in a 
more general context. In this paper, we construct self-similar 
solutions for "champagne fiows" with central voids in self- 
similar expansion. The inclusion of a central void not only 
makes our model more realistic, but also allows us to take 
into account stellar wind bubbles. 

We outline the model formulation of a general poly- 
tropic gas and present self-similar asymptotic solutions in 
section 2 and construct global solutions of "champagne 
fiows" in section 3. Section 4 provides solutions of self-similar 
"champagne fiows" with an expanding central void. In sec- 
tion 5, we discuss behaviours and astrophysical applications 
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of our novel solutions, and suggest other plausible forms of 
H II regions. Details of an invariant form of self-similar so- 
lutions in a conventional polytropic gas with the self-gravity 
ignored are summarized in Appendix A. 



By performing self-similar transformation ^ in equa- 
tions (HI — ® and introducing parameter q = 2{n -f 7 — 
2)/(3n — 2), we obtain two integral relations 



ax {nx ~ v) , 



(8) 



(9) 



2 SELF-SIMILAR POLYTROPIC FLOWS 

2.1 General Polytropic Formulation 

Dynamic evolution of a quasi-spherical general polytropic 
gas under self-gravity can be described by nonlinear hydro- 
dynamic equations in spherical polar coordinates (r, Q, 0), 



dp 1 d 2 . r^ 



dM dM 
dt dr 



dM 
dr 



( du du 



dp _ GMp 
dr r^ 



P = i^{r, t)p'< 



(1) 
(2) 
(3) 

(4) 
(5) 



where p(r, t) is the mass density, u(r, t) is the bulk gas 
radial flow velocity, M(r, t) is the enclosed mass within r at 
time t, p is the thermal pressure, G = 6.67 x 10~* dyne cm^ 
g~^ is the gravity constant. Equations H}, ([2]), ([3} describe 
the mass conservation, equation Q is the radial momentum 
equation and equation ([5]) is the general polytropic equation 
of state, in which 7 is the polytropic index and the coefficient 
K directly related to the 'specific entropy' depends on r and 
t. For a conventional polytropic gas, k is a global constant in 
space and time. More generally, we require the conservation 
of 'specific entropy' along streamlines, namely 



d_ d_ 

dt dr 



In — = 



(6) 



This s et of equations is the same as those of IWang fc Loul 
l|2008ll but without a completely random magnetic field. 

To reduce the nonlinear partial differential equations 
(PDEs) to ordinary differential equations (ODEs) for self- 
similar flows, we introduce the following transformation 



k^'H'^x 



p = 



kt 



u = k 

2n-4 



1/2 ,n-l 
' t V 



AttG 



M = 



ATvGf2 

fc3/2^3n-2^ 

(3n - 2)G 



(7) 



where a; is a dimensionless independent self-similar variable, 
A; is a dimensional parameter related to the polytropic sound 
speed making x dimensionless, v{x), a{x), /3{x), m{x) are di- 
mensionless reduced dependent variables of x only, and n is 
a key scaling index which controls the relation between r and 
X as well as various scalings of reduced dependent variables. 
We refer to v{x), a{x), f3{x), and m{x) as the reduced ra- 
dial flow speed, mass density, thermal pressure, and enclosed 
mass, respectively. Transformation ((Tjl is identical with that 
of Lou & Wang (2006). 



where there is no loss of generality to set the proportional 
coefficient (i.e., an integration consta nt) equal to unity i n 
integral ® for g / 2/3 or 7 7^ 4/3 (|Wang fc Loul I2OO8I ') . 
The special case of 7 = 4/3 corresponds to a relativistically 
hot gas as studied by Goldreich & Weber (1980) and Lou 
& Cao (2008). By setting 5 = 0, the general polytropic for- 
mulation reduces to the conventional polytropic case of a 
global constant k (e.g., Suto & Silk 1988; Lou & Gao 2006; 
Lou & Wang 2006; Lou, Jiang & Jin 2008) with n -I- 7 = 2. 
According to expression M{r, t) for the enclosed mass in 
transformation (O, we require 3n — 2 > and nx ~ v > to 
ensure a positive enclosed mass. The inequality 3n — 2 > 
will reappear later for a class of asymptotic solutions at large 
x. By equation we emphasize that for nx — v = at a, 
certain x* , the enclosed mass within x* becomes zero; we 
refer to this as a central void and x* is the independent sim- 
ilarity variable marking the void boundary which expands 
with time t in a self-similar manner. 

Combining all reduced equations above, we readily de- 
rive two coupled nonlinear ODEs for a' and v' as 

X{x,a,v)a' = A{x,a,v) , X{x,a,v)v' — V{x,a,v) , (10) 

where functionals X, A and V are defined by 

X{x, a,v)=[2-n+ (3n - 2)q/2] q1-"+3"9/2 

xx^''{nx — vy — {nx - v)'^ , 



A{x,a,v) = 2^c 



+(nx — v) 



qa^-^+^^^/'^x^^inx - v)"-^ 



+qa^-''+^"''^^x'^''-\nx - v)''-^3nx ~ 2v) 



V{x, a, v) = 2^Q (2-n+^q) q,-"+3"9/2^29 



X (nx — vY — {nx — v) (n — l)v + ^^-r§<^ 
+qa'--"+^"''^^x^''^\nx - v)''~''{3nx - 2v) 



(11) 



For a conventional polytropic gas of constant k with n-l-7 = 
2, we simply set q — in equation (|lip to derive 

a' _ ('^ ~ 1)^ + ''"3n-2)' - 2(a^ - ~ '")/^ 



a 



a{nx — vY — 'yet'' 

(3n-2) 



^ (n — l)av{nx — v) + — 2'^cP (x — v)/x 



(13) 



a(nx — w)2 — "ia~< 

Up to this point, our basic self-similar hydrodynamic for- 
mulation is the sam e as that of Lou & Wang (2006) and of 
IWang fc Loul (|2008l ) without a random magnetic field. 

For energy conservation, we define the energy density e 
and the energy flux density J as follows, 

pu^ GM p i 



+ 2^ 



(14) 
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GM 



2 p 



(15) 



where e is the energy density, J' is the energy flux density 
and i is the degree of freedom of an individual gas particle. 
The three terms in expressions (I14|l and (I15|) correspond to 
densities of the kinetic energy, the gravitational energy and 
the internal energy, respectively. With equations ((ij — (O 
and a globally constant k, we derive 



de 
dt 



+ 



J.2 Qj. 



dp 

'■IT 
or 



:(7-l) 



(16) 



for energy conservation, where V represents the net energy 
input. If the gas expands adiabatically or 7 = (i + 2)/i, 
then P = 0. Whether the gas locally gains or loses energy 
depends not only on the difference between 7 and (i + 2)/i, 
but also on the signs of dp/dr and u. 

2.2 Self-Similar Solutions 

An exact globally static solution, known as the singular poly- 
tropic sphere (SPS) takes the following form of 



,2-9 



-l/(n-3nij/2) 



"^i 2(2-n)(3n-2) J ' ^^^^ 

This is a straightforward generalization of the singular 
isothermal sphere (SIS; e.g., Shu 1977) and of the SPS for a 
conventional polytropic gas with g = (Lou & Wang 2006; 
Lou & Gao 2006). For a general SPS here, the mass density 
profile scales as p (x r"^/", independent of q parameter. 
For large x, the asymptotic flow behaviour is 

+ 2(2 - n)n^-iAi-"+^"^/^ 



(3n - 2) 



1-2/n 



-l/n 



+ ■ 



a = Ax 



-2/n 



+ ■ 



(18) 

(19) 



where A > Q and B are two constant parameters. As a; ^ 
+00, the mass density proflle scales the same way as SPS 
H17p above, independent of q parameter to the leading order. 
Since x +00 means f ^ 0"*", the asymptotic solution at 
large x is thus equivalent to initial conditions of the system 
at a finite r. Coefficients A and B are referred to as the 
mass and velocity parameters, respectively. A global solution 
with asymptotic behaviour (|18|l and (|19[) represents a fluid 
whose density profile scales initially similar to that of a SPS. 
Thus, by varying the scaling index n or the polytropic index 
7, we are able to model the initial density proflle with an 
index I = —2/n. Furthermore, as a physical requirement for 
plausible similarity solutions of our polytropic flow, both 
v{x) and a{x) should remain finite or tend to zero at large x. 
Hence, in cases of 2/3 < n < 1, corresponding to 2 < / < 3, 
A and B are fairly arbitrary, while in cases of 1 ^ n < 2 
(general polytropic), corresponding to 1 < I ^ 2, B should 
vanish for a finite radial flow velocity at large x. 

In the regime of a; — > 0^ , there exists a free-fall asymp- 
totic solutions for which the gravity is virtually the only 
force in action, and the radial velocity and the mass density 
profile both diverge in the limit of a; ^ 0"*". To the leading 
order, the free-fall asymptotic solution takes the form of 

1/2 



a{x) = 



(3n - 2)m(0) 



2x3 



(20) 



v{x) — 



2m (0) 



11/2 



(3n - 2)x 



(21) 



where constant m(0) represents an incre asing centra l point 
mass. Such solutions were first found bv IShrj ()l977h in the 
isothermal case, and were g eneralized to the conventional 
polyt ropic case (Cheng 1978: ILou fc Gaoll2006l : lLou fc Wan"3 
l2006f ) and to th e general polytropi c case with a random 
magnetic field bv lWang fc Loul l|2008l ). The asymptotic form 
does not depend on q. The validity of such solutions requires 
n > 2/3 and 7 < 5/3; the last inequality appears as a result 
of comparing various terms in series expansions. 

Another exact global solution, known as the Einstein- 
de Sitter (EdS) solution, exists in two cases of g = and 
g = 2/3. The EdS solution for a conventional polytropic gas 
of q = reads as 

2 2 2(n-2/3) 3 , , 

""3^' ""3' 3 ' ^22) 

(Lou & Wang 2006); for q = 2/3 and thus 7 = 4/3, this 
global EdS solution without and with a random magnetic 
field takes a slightly different form (Lou fc Cao 2008). With 
shocks, EdS solutions can be used to construct polytropic 
"champagne flows" with various upstream dynamics. 

From now on, we focus on the conventional polytropic 
case of g = with n-l-7 = 2. To simulate "champagne flows" 
in H II regions in star-forming clouds, we solve coupled non- 
linear ODEs (|12|l and (|13p subject to the inner boundary 
conditions at a; = 0, namely. 







(23) 



where ao is a constant. A series expansion yields the LP 
type of asymptotic solutions at small x in the form of 



(1-7) 




(24) 
(25) 



2 a 

V = —X 

3 157 

The isothermal counterp art of t his polytropic series solu- 
tion was obtained earlier l)Larsoni ri969a. b; Pension 1969a. ibl. 
Hunter 1977; Shu et al. 2002) with n = 1 and 7=1. Such 
LP type solutions may be utilized to construct champagne 
flows, if we ignore the central protostar as an approximation 
and assume the surrounding gas to be initially static (e.g., 
IShu et aLlliooi ). Physically, as a result of the gravity of the 
central protostar, the gas infall towards the very central re- 
gion may not stop, even when the protostar starts to shine 
and the 'champagne phase expansion' has occurred in the 
outer cloud envelope due to the photoionization and ultra- 
violet (UV) heating. Thus with central free falls (|20|l and 
pi[l as the downstream side of a shock, we can also con- 
struct global solutions for the dynamics of H II regions sur- 
rounding a nascent central massive protostar which involves 
free-fall materials. We can take LP type or EdS solutions as 
the downstream side of a shock and model classical cham- 
pagne flows for a conventional polytropic gas. We shall come 
to the possible scenario of inner free-fall solutions with an 
outer champagne flow. 

With the restriction n -f 7 = 2 for a conventional poly- 
tropic gas, the initial density profile with the scaling index 
n is directly linked to the polytropic index 7, depending on 
the energy exchange process in the gas (see eg I16|) . For a 
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general polytropic gas with g 7^ in contrast, the scaling 
index n and polytropic index 7 can be independently spec- 
ified; in particular, we can have 7 > 1 and n ^ 1 (this is 
impossible for a conventional polytropic gas). On the other 
hand, the initial mass density profile is affected by the star 
formation or other energetic processes before t = 0. Hence 
for I — — 2/n, we postulate that the energy exchange process 
remains largely unchanged after the protostar formation at 
t — 0. This is plausible because the ionization front travels 
relatively fast to large distances in a cloud during the "cham- 
pagne fiow" phase. Our general polytropic model allows in- 
equality 2/3 < n < 2, corresponding to 1 < / < 3, which 
covers so far the entire range of initial mass density profiles 
of H II regions. This range of mass density distribution has 
been obtained from radio obs ervations of cloud fragments 
and isolated dark clouds (e.g., lArguilla &: Goldsmith|[l985l : 
lMverslll985l ). 

Franco et al. (1990) reveals that the initial mass density 
profile index 3/2 < / < 3 (i.e., 2/3 < n < 4/3 in a poly- 
tropic model) leads to "champagne fiows" in clouds with 
weak shocks; and I > 3 corresponds to a "champagne flow" 
with strong and accelerating shocks. We then require pa- 
rameter n within the range 2/3 < n < 4/3, as we assume 
the cloud is 'density bounded'. With such range of parame- 
ter, we provide solutions of "champagne flow" with shocks. 
We emphasize in particular that there is one more degree of 
freedom to choose the velocity parameter B in asymptotic 
solution (|18|) and (|19|) in polytropic cases with 2/3 < n < 1 
(i.e. 2 < / < 3) than in the isothermal case. This leads to ma- 
jor differences of the polytropic champagne flow solutions we 
find as compared with those isoth ermal solutions of S hu et 
al. (2002) and Tsai & Hsu (1995). iFranco et all l|2000i ') have 
recently argued from the radio continuum spectra that UC 
H II regions have initial density gradients with 2 ^ / ^ 3; so 
we would consider primarily the n range of 2/3 < n < 1 or 
2 < Z < 3. In summary, molecular clouds with 3/2 < Z < 3 
have "champagne fiows" in self-similar manner, while those 
clouds with I ^ 3 have "champagne fiows" without simi- 
larity. For clouds with 2 < Z < 3, there is more than one 
parameter to specify initial dynamic flows. 

Finally, we need to include shocks at proper places in 
LP type solutions, EdS solutions or free-fall solutions to 
match with appropriate asymptotic solutions at large x to 
determine relevant coefficients. The shock jump conditions 
between the downstream and upstream variables are de- 
termined by the mass conservation, the radial momentum 
conservation and the energy conservation. With these three 
equations for shock conditions, one can determine upstream 
self-similar variables {xsu, Qu, Vu) uniquely from down- 
stream self-similar variables {xsd, Od, Vd) or vice versa. 
The two subscripts d and u here denote the downstream 
and upstream variables, respectively. Detailed formulation 
and procedure of self-similar shocks can be found in Section 
5 of Lou & Wang (2006). All solutions in this paper are ob- 
tained by solving coupled nonlinear ODEs p2p and (|13p for 
a conventional polytropic gas with n + 7 = 2. 



3 POLYTROPIC CHAMPAGNE FLOWS 

In cases of 71 < 1, there is a range of shock positions (or 
speeds) for a specified downstream solution with a fixed den- 



sity at the centre ao at a: = 0, corresponding to different 
asymptotic fiow behaviours at large x on the upstream side. 
We now discuss two situations: first, for cases with a fixed 
value of ao, we adjust shock positions for a specified LP 
type solution and observe the relation between shock posi- 
tions and asymptotic flow behaviours at large x. Secondly, 
for cases with a fixed shock position, we alter the value of qo 
and examine the variation of upstream conditions. It is ex- 
pected to set certain limits on relevant parameters for poly- 
tropic "champagne flows" to exist. As a series of examples, 
we choose the scaling index n = 0.9. Numerical explorations 
have also been performed for cases of n = 0.7 and n = 0.8 
and the results are qualitatively similar. 

3.1 Cases with a Fixed ao Value 

With a flxed value of ao for x ^ 0"*", one can uniquely 
determine a LP type solution by a standard numerical in- 
tegration. Such a LP type solution will encounter the sonic 
critical curve at a certain ajmax uniquely corresponding to 
ao. It is natural to consider possible hydrodynamic shock 
positions Xsd < a^max on the downstream side of shock front. 
One such example of n = 0.9 and ao = 1 is shown in Fig- 
ure [1] with relevant parameters summarized in Table [T] To 
model "champagne fiows" with outward expansions at large 
radii in molecular clouds, we provide the following analysis. 
In principle, there are two other conditions giving rise to two 
minima of x^d in order to obtain an outfiow in the upstream 
side of a shock. First, our extensive numerical explorations 
reveal that there exists a a;mini and for Xsd < a^mini, the 
upstream shock position Xsu becomes complex by the shock 
conditions. This only happens when we attempt to obtain 
upstream variables from downstream variables. For a real 
Xsu, downstream variables {xsd, ad, Wd) should satisfy 



(l-7)a]-i+2(nx,d-«d)'>0 



(26) 



For 7 < 1 (unphysical) and 7 = 1, inequality (|26p is readily 
satisfied; for 7 > 1 or n < 1, this condition does not always 
hold. Algebraic manipulations give the downstream Mach 
number in the shock reference framework A^d as 



^2 ^ {nxsd - Vd) 



7-1 



(27) 



Inequality (|26|l imposed on a subsonic downstream Mach 
number is 1 > Aid > (7 ~ l)/(27). The downstream Mach 
number and the upstream Mach number is related by 



Ml- 2 + (7-l)-M^ 



27.M g-(7-l) • ^^^^ 

This relation was provided bv lLou fc Cad (|2008l ) for a rela- 
tivistically hot gas and is proven valid in our model consid- 
eration. With the possible range of upstream Mach number 
1 < Mu < +00, we then have the limit on Md shown above. 
As we integrate LP or EdS solutions from a; = with a = ao 
and V = 0, solutions do not satisfy inequality l|26p when x 
remains sufficiently small and Xmini is the minimum value of 
X satisfying (|26|l . This value of a^mini is uniquely determined 
by the ao value. Therefore, for a fixed LP type solution 
around small x, polytropic shocks can be constructed with 
a downstream shock in the range of a;inini < Xsd < Xm^x, 
and across such a shock, the LP type solution at small x 
can be matched with different asymptotic fiows at large x. 
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Systematic numerical explorations reveal that the up- 
stream velocity increases monotonically with the increase of 
3^sd) as shown by the variation trend of the B parameter 
(see Table [l]). There is thus another critical value imposed 
on Xsd, denoted by a;min2- For Xsd > a;min2, the upstream 
solution matches to an asymptotic solution at large x in the 
form of (|18|) and (|19|) with B > Q, referred to as an outflow. 
Complementarily with Xsd < Xmi-n2, the upstream solution 
matches to an asymptotic solution with B < 0, referred to as 
an inflow. As B varies continuously and monotonically with 
XaA, for Xsd = Xyaiyi2 the upstream solution corresponds to 
an asymptotic solution with _B = 0, which describes a breeze 
or a contraction in association with "champagne flows" . Ac- 
cording to asymptotic expressions (|18p and (|19|) with g = 0, 
the breeze or contraction correspond to slow outward or in- 
ward flows. To obtain a breeze, we need a mass parameter 

For the specific case oi A = As, either a breeze or a con- 
traction reduces to SPS solution (I17p . With n < 1, there are 
three possibilities in general. First, Xmini > Kniin2 for the 
allowed range to construct shocks, the upstream solutions 
always correspond to outfiows. Secondly, Xmmi < a;min2 for 
the allowed range to construct shocks, it is possible to ob- 
tain outflows, inflows and breezes or contractions for the up- 
stream side. Thirdly, if any of a;mini or a;min2 exceeds a;inax, 
a global champagne flow is not allowed. For the isother- 
mal case of n = 1, we may set B equal to zero; asymptotic 
breezes or contractions on the upstream side are allowed, 
and the only al lowed value of the downstream shock position 
IS Xsd — ^min2. IShu et al.l (|2002l ) indicated that the isother- 
mal shock position is uniquely determined by the value of 
Qo, which is consistent with our mo re general analysis here. 
The unique shock position found bv lShu et al.l (2002) corre- 
sponds to the a;min2 above. 

The conventional scenario for "champagne flows" would 
require the entire fluid to expand outw ard. Numerical sim- 
ulations on "champagne flows" (e.g., iTenorio-Tagle et al.l 
[1986 ) assume that at t = 0^ the central star is formed and 
the surrounding cloud is initially at rest. For t > 0'^ , the 
fluid is photoionized and heated by the ultraviolet radia- 
tion from the central star and expands. In this scenario, we 
would require an expanding upstream flow in order to model 
"champagne flows" . However, since solutions with an asymp- 
totic inflow or contraction as upstream part may also exist, 
the outer part of H II regions can also have inward veloci- 
ties. In fact, during the star formation even some time after 
the star formation, the surroun ding cloud may cont inue to 
collapse towards the centre (e.g.. iFatuzzo et al.ll2004h . With 
the core nuclear burning of the central protostar, the sur- 
rounding gas is ionized and heated, and the inner part of the 
fluid starts to expand, while the outer part continues to fall 
inwards. Those solutions of the LP type as the downstream 
side of a shock and an asymptotic inflow or contraction for 
the upstream side correspond to this scenario just described; 
and such global solutions are referred to as Inner Shock Ex- 
pansions in a Collapsing Envelope (ISECE). 

For the case of n = 0.9 and ao = 1, we have deter- 
mined Xmini = 0.95, 2;min2 = 2.287 and Xmax = 3; the shock 
range 2.287 < < 3 gives sensible classical poly tropic 
"champagne flow" solutions and 0.95 < Xsd < 2.287 gives 
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Figure 1. The reduced mass density a(x) (top) and the reduced 
radial flow velocity v{x) (bottom) for global "champagne flow" 
solutions in cases with n = 0.9 (thus 7 = 1.1) and ao = 1- In 
both panels the dashed curve represents the sonic critical curve 
and in the bottom panel the dotted line is i; = 0. The downstream 
solution is connected with various upstream solutions with light 
solid curves. The downstream solution is integrated numerically 
from — » 0"*" with a LP type solution. In both panels, the up- 
stream solutions from top to bottom correspond to Xsd=3, 2.923, 
2.773, 2.623, 2.473, 2.287, 2.2 and 2 respectively. We note that 
the upstream solution with Xsd = 2.287 corresponds to a breeze. 
Relevant parameters are summarized in Table [Tl 

the ISECE solutions as shown in Figure [T] For n — 0.9, we 
have As — 2.042; with Xsd = a;min2 and A = 1.536 < As, 
the asymptotic solution is a breeze and should be consid- 
ered as a classical champagne flow. Given other parame- 
ters the same, the situation oi A > As would give rise to 
an asymptotic contraction. The shock location and shock 
speed can be determined once Xsd is specified. The dimen- 
sionless shock position in the self-similar variable represents 
the shock strength and velocity in dimensional form. The 
shock velocity reads drs/dt = nky^Xsdt"~^ ■ The outgoing 
shock slightly decelerates (for n slightly less than 1) and the 
shock velocity is proportional to Xsd- 

As shown in Figure [T] and Table [TJ with increasing Xsd, 
the upstream variables at the shock front Vu and Qu increase, 
and the two parameters A and B of asymptotic solutions 
(|18|) and (|19p also increase. This is a fairly common feature, 
observed in all other cases that we have studied numerically. 
Different shock positions match with different asymptotic 
solutions at large x. Once A and B are specified, the shock 
position Xsd is uniquely determined. As A is the mass pa- 
rameter and B is the velocity parameter, Xsd, proportional 
to the shock velocity and strength, is determined not oifly 
by initial mass density but also by the initial motion. This 
differs from the isothermal case when B — with the mass 
parameter A determining the shock behaviour. We expect a 
faster and stronger shock with a higher initial speed. 

We further identify two sub-types of such ISECE solu- 
tions: (i) the upstream side has an outward velocity near the 
shock and an asymptotic inward velocity far from the centre 
(e.g., solution with Xsd ~ 2.2 in Figure [TJ; (ii) the upstream 
side has an inward velocity everywhere (e.g., solution with 
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Xsd = 2 in Figure [TJ . For the type (i) solutions, there is a 
stagnation point Xstg where the radial flow velocity vanishes; 
for the solution with Xsd = 2.2, Xstg ~ 3. With self-similar 
transformation ([7|, this stagnation point rstg = k^^'^t" Xatg 
travels outward with time in a self-similar manner. 

For ISECE solutions, we envision that such solutions 
correspond to the situation where a star starts to burn, ion- 
izing and heating the surrounding medium as the gas falling 
continues. The gas infall and collapse are indispensable in 
star formation. If the nascent star ionizes the whole resid- 
ual gas sufficiently fast, the outer gas may still possess an 
inward momentum. A champagne shock runs into the in- 
fall gas, deposits outward momentum and accelerates the 
outer gas. If a shock is sufficiently strong, we expect type 
(i) solutions, e.g., the gas immediately outside a shock flows 
outward, and the stagnation point travels outward propor- 
tional to the self-similar expansion of the shock. If a shock is 
sufficiently weak, we expect type (ii) solutions, e.g., the gas 
flows inward outside the shock front. This ISECE scenario 
is expecte d to occur in certain H II regions (e.g., Shen & 
Lou 2004: iFatuzzo et al.l (|2004l )). 



3.2 Cases of a Fixed Dimensionless Shock Position 

The variation of LP type solutions with different qq and the 
same a;sd is shown in Figure (2] For ao > 2/3, a decreases 
with increasing x while for ao < 2/3, a increases with in- 
creasing X. With a larger ao the LP type solution encounters 
the sonic critical curve earlier; therefore ao cannot be too 
large, otherwise the LP type solution encounters the sonic 
critical curve before reaching the preset shock position Xsd- 
One case of n = 0.9 and Xsd = 3 is shown in Figure [2] with 
parameters given in Table [S] Here, we have a^d < 2.5. 

As shown in Figure [2] and Table [51 with the increase of 
ao, the upstream conditions at the shock front and the 
mass parameter A of the asymptotic solution increases, but 
Hu and B do not have a steady trend. As ao increases from 
10~^ to 0.04, Vu and B decrease, while as ao increases from 
0.04 to 2.5, Vu and B increase (see Tabled . This shows that 
Vu and B are correlated. The minimum values of Vu and B 
are associated with ao — 0.04. Here we show that with a 
prefixed downstream shock position x^d, it is possible that 
neither breeze nor contraction solutions are allowed. With a 
Xsd as large as 3, all the asymptotic upstream sol utions cor- 
respon d to outflows. In the isothermal analysis of lShu et al.l 
l|2002h . the mass parameter A tends to zero with ao — > 0"*". 
In polytropic cases, we obtain similar results. According to 
series expansions (|24p and (|25|l . if setting ao = exactly, 
the integration gives a trivial solution of a = 0. 

Base d on numer i cal ex plorations for ao ^ 0^ in isother- 
mal cases. Ishu et al.l (|2002| ) argued that the self-gravity may 
be neglected for cases with small central density, and devel- 
oped another self-similar transformation, viz., the so-called 
invariant form, in order to model the initial mass density 
profile other than that of a SIS (i.e., p oc r~^). The initial 
mass density profile p oc r~' , where index I does not nec- 
essarily equal to 2, can be described by the invariant form 
when the self-gravity is ignored. We perform a similar re- 
duction with the invariant self-similar transformation for a 
conventional polytropic gas without the self-gravity (see Ap- 
pendix A), and show that with n 7^ 1 (i.e., non-isothermal 
cases), the power index / must be equal to 2/n for a self- 



n=0.9, X =3 

sd 



S 2 





2.5 
2 

1.5 
1 

0.5 
























/ a =2.5 a=OM a =0.5 a=^0'^ ^a=W^ 







10 



Figure 2. The reduced mass density ol(x) (top) and the reduced 
radial flow velocity v(x) (bottom) for global polytropic "cham- 
pagne flow" solutions in cases n = 0.9 (thus 7 = 1.1) and Xgd = 3. 
In the top panel, the horizontal dotted line stands for q = 0, sep- 
arating the top panel into two parts; the vertical scales in these 
parts are different. In the upper part a{x) is presented linearly, 
while in the lower part log[«(a;)] is presented. In the bottom panel, 
the dotted line is the sonic critical curve. Shocks at a; ~ 3 as 
the discontinuity in solutions are shown. The solid curves on the 
downstream side of shocks are LP type solutions with ao=0.5, 
1, 1.5, 2 and 2.5 (from bottom to top in the top panel and from 
top to bottom in the bottom panel). The solid curves on the 
right show respectively the corresponding upstream solutions ap- 
proaching different asymptotic solutions at large x (from bottom 
to top in the top panel and from bottom to top in the bottom 
panel). The dashed curves on the downstream side of shocks are 
LP type solutions with ao = 10"^, 10"*, 10"^, and 10"^ (from 
bottom to top in the top panel and from top to bottom in the bot- 
tom panel). The dashed curves on the right show respectively the 
corresponding upstream solutions approaching different asymp- 
totic solutions at large x (from bottom to top in the top panel and 
from top to bottom in the bottom panel). The thick solid curve in 
both panels represents the solution with ao = 0.04, which has the 
lowest upstream speed and the minimum B at large x. Relevant 
parameters are summarized in Table |2] 
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Table 1. Data Parameters of Global Polytropic "Champagne Flow" Solutions in Cases with n = 0.9 and ao = 1 



A 


B 














0.9902 


-0.7699 


2 


0.8416 


1.3019 


2.0693 


0.2005 


-0.3012 


1.3501 


-0.2267 


2.2 


0.8172 


1.4266 


2.2446 


0.2404 


0.1004 


1.5357 





2.2869 


0.8068 


1.4805 


2.3234 


0.2586 


0.2599 


2.0109 


0.4767 


2.4732 


0.7848 


1.5956 


2.4964 


0.2993 


0.5789 


2.474 


0.8577 


2.6232 


0.7678 


1.688 


2.6387 


0.3336 


0.8172 


3.0224 


1.2416 


2.7731 


0.7516 


1.7803 


2.7832 


0.369 


1.0425 


3.6696 


1.6323 


2.9231 


0.7363 


1.8728 


2.9292 


0.4054 


1.2568 


4.0455 


1.8367 


3 


0.729 


1.9204 


3.0047 


0.4243 


1.3629 



similar form. This mass density profile with a scaling index 
I = 2/n\s the same as that in a SPS and in asymptotic solu- 
tions at large x. In summary, the freedom of choosing I in the 
invariant form disappears for non-isothermal cases. From an- 
other perspective, since the index n ranges in 2/3 < n < 4/3, 
self-similar polytropic champagne flows can model the ini- 
tial mass density profile with 3/2 < Z < 3. In other words, 
the objective to model the initial density profile other than 
Z = 2 is naturally fulfilled, without the necessity of drop- 
ping the self-gravity. The clear advantage of our polytropic 
approach is that the self-gravity is included in the model 
consideration. Therefore, to apply our polytropic solutions, 
ao — > 0^ is no longer required. From transformation (0, pa- 
rameter ao is tightly linked with the central mass density po 
and timescale t, and should have different values in different 
situations. Therefore polytropic champagne fiow solutions 
are adaptable to a much wider range of astrophysical cloud 
systems. Moreover, the ao 0"*" cases in our polytropic 
framework can be approximated by a central "void" as dis- 
cussed in the following section. With a central "void" , we are 
able to neglect the gravity of the central region where the 
density is sufficiently low and still consider the self-gravity 
of the outer more dense gas medium. 

In the isothermal case of lShu et al.l(|2002l ). the solution 
in which the outer part is a static SIS represents a limiting 
solution defining the maximum value of ao. In the polytropic 
cases with n = 0.9, we can also identify such a limit by re- 
quiring B = Q for the upstream asymptotic solutions, such 
that for a preset ao the downstream shock position Xsd is 
uniquely determined. A family of such "champagne flow" so- 
lutions with asymptotic upstream breezes or contractions is 
shown by solid curves in Figure |3] with parameters summa- 
rized in TableO With a gradual increase of ao, the upstream 
side varies from outward breezes, to a SPS, and to an inward 
contraction; this trend leads to a maximum ao if we define 
an outward breeze for classical "champagne fiows". Here, the 
critical value ao = 3.13 fo r an upstream SP S corresponds to 
the upstream SIS limit in lShu et all (|2002| ) (i.e., ao = 7.9 in 
the isothermal case). Naturally, this critical value depends 
on the choice of n. For ao > 3.13, an upstream asymptotic 
contraction or an ISECE solution appears (dashed curve in 
Fig. |3}. Thus, isothermal champagne flows form a special 
family with B = and n = 1, referred to as breeze cham- 
pagne fiows. With 2/3<n<lor2<Z<3, there are 
many more physically possible champagne flows, for which 
the upstream side corresponds to asymptotic outflows at 
large x. The physical meaning of this special solution whose 
upstream side is the outer part of a static SPS is clear: the 



outer envelope of gas is initially in a hydrostatic equilibrium, 
and the expanding shock created by the UV photoionization 
travels into the static envelope; on the downstream side of 
this expanding shock, the gas is heated to high tempera- 
tures. As the static SPS relies on the scaling parameter n or 
the polytropic index 7, we expect one single solution for a 
preset n. From Table Owe see clearly that with the increase 
of ao, the downstream shock position to obtain a breeze in 
the upstream, or a;min2 decreases. Meanwhile, numerical ex- 
plorations suggest that a;mini increases with increasing ao. 
Hence we expect for a sufficiently large ao, a;min2 would be- 
come less than a;mini to forbid ISECE solutions. 



4 SIMILARITY POLYTROPIC "CHAMPAGNE 
FLOWS" WITH CENTRAL VOIDS 

We now establish and analyze a new class of "champagne 
flow" solutions with voids surrounding the centre. We ex- 
tend solutions from x = to the dimensionless self-similar 
expanding boundary 2;* of a void, inside of which there is 
no mass, i.e., m* = m{x*) = 0. In our notations, super- 
script * attached to variables indicates variables on the void 
boundary x* . By expression ([S]), we have v* — nx* . The void 
boundary conditions are 

a = a* , V — nx* , at x = x* . (30) 

A Taylor series expansion to the first order around the void 
boundary a; = x* of ODEs and yields 

v{x) ^ nx* + 2{1 - n){x - X*) ^ , (31) 

a{x) = a* + "(l-'^) ta*)''x*(x - x*) + ■ ■ ■ . (32) 

7 

Series expansions (|31|l and (I32|l are conspicuously different 
from series expansions (|24|l and (|25|l . The proportional co- 
efficient of the asymptotic velocity is no longer 2/3 but de- 
pends on the scaling index n. By translating the origin, the 
dynamic flow behaves differently. For n < 2, solutions given 
by expression (|31[) are locally above the line nx — v = 0, 
indicating a positive enclosed mass. Numerical integrations 
reveal that the solutions are always above the line nx — v = 
thereafter. Here, we model a gas flow with a central void in 
the presence of self-gravity and thermal pressure, directly 
relevant to galactic subsystems such as H II regions. 

One can readily obtain the downstream portion of 
"champagne flow" solutions by numerically integrating cou- 
pled nonlinear ODEs (|12|) and ((13} from the void boundary 
X* with asymptotic expansions l|3ip and (|32|) . All such void 
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Table 2. Data Parameters of Global Polytropic "Champagne Flow" Solutions in Cases with n = 0.9 and x^^ = 3. 
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0.00001 


0.0002532 


1.145 


3 


0.0001 


2.3469 


3.0609 


0.000030184 


1.5294 


0.0001 


0.001735 


0.9499 


3 


0.0007 


2.3152 


3.064 


0.0001993 


1.3959 


0.001 


0.01243 


0.7205 


3 


0.0049 


2.2821 


3.0693 


0.001372 


1.2353 


0.01 


0.0912 


0.4713 


3 


0.0358 


2.2449 


3.0746 


0.009752 


1.0536 


0.04 


0.3045 


0.3899 


3 


0.1154 


2.2135 


3.0715 


0.03192 


0.9641 


0.5 


2.4053 


1.1453 


3 


0.6048 


2.0444 


3.0197 


0.2556 


1.1567 


1 


4.0455 


1.8367 


3 


0.729 


1.9204 


3.0047 


0.4244 


1.3629 


1.5 


5.5009 


2.406 


3 


0.7512 


1.8166 


3.008 


0.5556 


1.5059 


2 


6.9587 


2.9308 


3 


0.7414 


1.7274 


3 


0.6606 


1.6083 


2.5 


8.5684 


3.4633 


3 


0.7588 


1.562 


2.8444 


0.7086 


1.4913 



Table 3. Parameters of Global Polytropic "Champagne Flow" Solutions with n = 0.9 and Upstream Breeze or Contraction {B = 0). 



0.2209 


0.7550 





2.6250 


0.3571 


1.8473 


2.6873 


0.0990 


0.5153 


2/3 


1.3240 





2.3986 


0.6667 


1.5991 


2.4419 


0.2045 


0.3407 


1 


1.5357 





2.2869 


0.8068 


1.4805 


2.3234 


0.2586 


0.2599 


2 


1.8643 





2.0733 


1.0783 


1.2563 


2.0996 


0.3734 


0.1069 


3.13 


2.0420 





1.9264 


1.2773 


1.1036 


1.9475 


0.4643 





5 


2.1924 





1.7706 


1.5091 


0.9429 


1.7873 


0.5746 


-0.1163 



solutions encounter the sonic critical curve, and the lower 
the value of q* is, the later (or at larger x) the void solution 
encounters the sonic critical curve. In order to match with 
asymptotic solutions of finite density and velocity at large x, 
void solutions must either cross the critical curve smoothly 
or connect to another branch of solutions via shocks. To 
model "champagne flows", we need to construct shocks to 
obtain global solutions. Global void solutions that cross the 
sonic critical curve smoothly will be discussed in Lou & Hu 
(2008, in preparation). We now present the solutions for 
large and small a*, respectively. A family of semi-complete 
"champagne flow" solutions with n = 0.9, void boundary 
X* = 1 and a* = 5 is constructed by varying the self-similar 
shock position as shown in Figure [l] Complementarily, an- 
other family of solutions with a* = 10"* is also constructed 
and shown in Figure [5] Both outflow, inflow and breeze or 
contraction as the upstream side are presented. Relevant 
parameters of these void "champagne flow" solutions are 
summarized in Table |4] 

Similar to no void cases, a decreases with increasing x 
for large a*, so the density on the void boundary is a local 
maximum (see Fig. |4]); while for small a*, a increases with 
increasing x, so the density maximum is at the downstream 
side of the shock. The latter corresponds to a shell-like struc- 
ture in self-similar expansion (see Figure [Sj. 

With downstream void solutions and upstream outflow 
and breeze solutions connected by shocks, we establish semi- 
complete polytropic "champagne flow" solutions with cen- 
tral voids. Similar to LP type solutions with shocks, there 
are also one maximum Xmax and two minimum limits Xmini 
and Xniin2 imposed on the downstream shock position a^^d 
in order to obtain "champagne flow" solutions. System- 
atic numerical explorations for cases of n = 0.7, n — 0.8 
and n = 0.9 show that in general a;min2 > a;mini. For 



2;mini < X sd < a;miii2, a void solution can be matched with 
an asymptotic inflow to produce ISECE solutions, while for 
Xsd > a;min2, a central void solution can be matched with an 
asymptotic outflow to produce "champagne flow" solutions. 
For Xsd = a;niin2, the upstream corresponds to a breeze or a 
contraction with B = 0. The analysis here parallels that in 
the previous section for cases without central voids; in par- 
ticular, the parameters a^mini and a;min2 are determined not 
only by n and a*, but also by the expanding void bound- 
ary X* . Numerical explorations suggest that for a certain n, 
with the increase of a*, a^mini increases and a;niin2 decreases. 
Hence for a sufficiently large a* , we expect Xniin2 < ajmini 
for which ISECE solutions are not allowed. This is consistent 
with polytropic cases without central voids. 

5 ANALYSIS AND DISCUSSION 

5.1 Comparison with Numerical Simulations 

To adapt our self-similar solutions for modelling an astro- 
physical cloud system, we need to first specify the parameter 
k related to the sound speed squared. By varying k, one can 
model clouds of different scales using a single self-similar so- 
lution. Parameter k is determined by the thermodynamic pa- 
rameters, including thermal pressure p, mass density p and 
temperature T. A useful relation for a conventional poly- 
tropic gas derived from transformation ((7} is 



where /i is the mean molecular mass of gas particles. 
For a typical value provided by the classification of 
iHabing fc Israel (|1979D . the UC H II regions have an electron 
number density Tie > 3000 cm~^ (corresponding to a mass 
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Table 4. Polytropic "Champagne Flow" Solutions with a Central Void inside x* = 1 in the Case of n = 0.9 



a* 




A 


B 






'"d 








10" 


-4 


0.000106 


-0.3813 


1.8 


1.4734 X 10""' 


1.3322 


1.9058 


2.7895 X 10"^ 


0.1058 


10" 


-4 


0.000191 





2.068 


1.8377 X 10-'' 


1.5386 


2.1442 


4.2673 X 10"^ 


0.4893 


10" 


-4 


0.0014 


1.0833 


3.088 


0.000517 


2.398 


3.152 


0.00015 


1.501 


10" 


-4 


0.0991 


3.5315 


5.701 


0.0114 


4.701 


5.839 


0.0032 


3.669 


5 




2.2512 


-1.0512 


1.7 


2.8906 


1.0956 


1.8424 


0.4646 


-1.2708 


5 




2.9388 





1.818 


2.4535 


1.1047 


1.8777 


0.5978 


-0.5632 


5 




5.2894 
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Figure 3. "Champagne flow" solution with a LP type down- 
stream and the upstream part as a breeze or contraction for 
n = 0.9 (i.e., 7 = 1.1). The top panel shows the density and 
the bottom panel shows the velocity. In both panels the dotted 
curve is the sonic critical curve. The solutions are integrated from 
ao = 0.2209, 2/3, 1, 2, 3.13 and 5 (from bottom to top in the top 
panel and from top to bottom in the bottom panel), and the 
downstream shock positions Xj^j are carefully chosen to let the 
upstream solution correspond to the asymptotic solutions with 
B = 0. Relevant parameters are summarized in Table |3] For the 
solution with 00 = 3.13, the corresponding upstream is a SPS. 
For ao > 3.13, the upstream contracts (dashed curve), while for 
00 < 3.13, the upstream is a breeze. For a breeze "champagne 
flow", we require ao ^ 3.13 (solid curves). 



density p>5xl0- g cm- for a fully ionized hydrogen 
gas), and the Compact H II regions have 1000 < < 3000 
cm~^ (corresponding to 1.7 x 10^^^ < p < 5 x 10"^^ g cm"^ 
for mass density p). Typically, the temperature of H II re- 
gions is of order of ~ 10* K. For a fully ionized hydrogen gas, 
we assume that fi = mp/2 where nip is the proton mass. The 





Figure 4. The reduced mass density a{x) (top) and the reduced 
radial flow velocity v{x) (bottom) for semi-complete "champagne 
flow" solutions with a central void inside x* = 1 in the case of 
n = 0.9 and density on the void boundary a* = 5. In both panels, 
the dashed curve stands for the sonic critical curve. In the bottom 
panel the dotted line stands for the line 1; = 0. The solid curve 
on the upper left of the sonic critical curve is the downstream 
void solution and the solid curves on the lower right of the sonic 
critical curve are the corresponding upstream solutions with the 
downstream shock position x^,^ = 1.7 (inflow), 1.818 (contraction) 
and 2 (outflow). In this case, Xmin2 ~ 1.818 and 3^minl 

= 1.408. 

Numerical data about these solutions are tabulated in Table |4] 



value of the polytropic index 7 does influence very much the 
resulting k and thus k should be evaluated specifically. For 
nearly isothermal cases, with relation (|33|l . we estimate k 
for UC and Compact H II regions to be fc ~ lO'^'^ ~ 10^^ cgs 
unit. One should be aware that k and thus k vary with the 
gas temperature and density in a cloud. Here, we presume 
a constant k to convert self-similar variables to real space 
variables as a first approximation. 

We now compare our self-similar solutions of quasi- 
spherical symmetry with previous numerical simulations. 
iTenorio-Tagle et al.l l|l986l) performed a numerical study for 
a similar scenario as ours, i.e. a nascent central massive pro- 
tostar ionizes and heats the ambient neutral gas and then 
leads to "champagne fiows". In their simulation, the radia- 
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Figure 5. The reduced mass density a(x) (top) and the reduced 
radial velocity (bottom) for semi-complete "champagne flow" 
solutions with a central void inside x* = f in cases with n = 0.9 
(thus 7 = 1.1) and a density on the void boundary a* = 10~*. In 
both panels, the dashed curve represents the sonic critical curve. 
In the bottom panel, the dotted line stands for the line v = 0. The 
solid curve on the left side is the downstream void solution and 
the solid curves on the right side are the corresponding upstream 
outflow solutions with the downstream shock position x^d = 1.8 
(inflow), 2.068 (breeze), 3.088 (outflow) and 5.701 (outflow). Here, 
= 2.068 is the limit to ensure an asymptotic outflow, i.e., 
3;min2 ~ 2.068. In this case, x^i^i = 1.261. Numerical data for 
these solutions are tabulated in Table |4] 



tive cooling rate is assumed to be low and thus our poly- 
tropic approach may be applicable. Franco et al. (1990) gave 
an analytical model for the formation and expansion of H II 
regions and have their solutions compared with simulations. 
We intend to demonstrate that our self-similar polytropic 
analysis of the problem gives qualitatively similar results. 

In the simulation o f iTenorio-Tagle et al.l (|l986D . compu- 
tations were carried out following the progressive ionization 
of a diffuse gas and subsequent dynamical evolution, in a 
globular cluster soon after the star formation has been initi- 
ated. The residual gas is initially in a hydrostatic equilibrium 
in the gravitational field at a uniform temperature and all 
the ionizing UV radiation comes from stars at the very cen- 
tre of the cluster. We emphasize that only the cases in which 
the gas is fully ionized correspond to our conventional "poly- 
tropic champagne flow" model. In the simulation, the initial 
mass density at the centre po, initial temperature To, and the 
UV photon flux F from the central star completely deter- 
mine the evolution in spherical symmetry. To translate these 
parameters into our self-similar form of solutions, we first ex- 
plore the time evolution of a "champagne flow" shock. In our 
self-similar model, the shock radius obeys rs — k^^^Xst". 
With such a self-similar formula, we can fit the scaling index 
n and k^^'^Xs to shock positions obtained by the n umerical 
simulation for case F of lTenorio-Tagle et al.l (|l986D . and the 
result comparison is shown in Figure |6] with relevant param- 
eters in caption. We see an almost perfect fit, suggesting that 




5.5 

Log time [yrs] 

Figure 6. Shock position evolution with time 4 in a "cham- 
pagne flow" for our self-similar model (solid li ne) and the nu- 
merical simulation of iTenorio-Taele et al.l (Il986f) (asterisks). For 
the simulation, the adopted parameters are central initial den- 
sity Po = 2 X 10"-^^ g cm~3, initial temperature Tg = 3000 K, 
and the ionization UV flux F = 2 X 10^^ photons s^^. For the 
self-similar model, the best fitted parameters are n = 1.0583 and 
log(A;^/^a;s) = 6.0356 in cgs unit. 



the dynamical evolution of "champagne flows" approaches 
a self-similar form. This fitting gives a value of n = 1.0583 
and log{k^^^Xs) — 6.0356 in cgs unit (strictly speaking, we 
need a general polytropic gas for 7 > 1). The value of n is 
fairly close to the isothermal case of n = 1, consistent with 
the model analysis of the num erical simulation that the gas 
evolution is nearly isothermal (iTenorio-Tagle et al.lll986l ). 

We now generate a global "champagne flow" solution 
grossly comparable to figure 2 of lTenorio-Tagle et al.l (|l986l ) 
by fitting parameters. We first choose a central reduced den- 
sity ao '~ 1 X 10~^ such that the initial central density yields 
the value used in the simulation. We expediently choose 
n = 0.9 (hence the parameter 7 = 2 — n = 1.1) to model the 
initial density profile I = —2/n — —2.22. Note that parame- 
ter n for the "champagne flow" solution is slightly different 
from the value we obtain from the fitting. With po and To 
specified in the simulation, we estimate kd = 3.6 x 10^^ 
cgs unit with expression p3p . We still have the freedom 
to require the shock traveling to Ts = 2.51 x 10^^ cm at 
t — 1.3 X 10^ yr, giving Xsd = 1.86. With these parameters, 
we model "champagne flows" in diffuse H II regions with 
radius r up to lO'^^ cm (i.e., ~ 300 pc). The full solution is 
shown in Figure [T] The timescale of 1.3 x 10^ yr is regarded 
as the duration of the formation phase and the initial time 
for a "champagne flow", and the timescale 5.1 x 10^ yr is 
regarded as the lifetime of a "champagne flow" . 

The orders of magnitude of all variables are consistent 
with typical values; e.g., the expansion velocity is several 
tens km s~^ and the temperature is about ~ 10* K. The 
enclosed mass at r = 10^^ cm (~ 300 pc) is about 850Mq, 
consistent with the value of ~ 800A/q given by the numerical 
simulation and with the typical value for diffuse H II regions. 
The enclosed mass does not vary with time t, conflrming the 
cut-off radius chosen at r = lO^'^ cm. As time evolves, the 
central number density decreases from 10~°'^ to 10~^'^ cm~^ 
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Figure 7. Self-similar "champagne flow" solutions for radius up 
to 10^^ cm {~ 300 pc) at time t = 1.3 X 10^ yr (solid curves) and 
t = 5.1 X 10^ yr (dashed curves). From top to bottom, the pan- 
els show number density, velocity, pressure, enclosed mass and 
temperature of the gas, respectively. The self-similar shock so- 
lution is obtained with n = 0.9, 7 = 1.1, ao = 1 X 10~^ and 
downstream shock position x^^ = 1.86. The downstream sound 
scaling factor is 3.6 X 10^^ cgs unit, and the upstream sound 
scaling factor ku is 3.38 X 10^^ cgs unit. The self-similar vari- 
ables on the downstream side of the shock are {xg^, a^, v^) = 
(1.86, 2.80 X 10~^, 1.37), and the corresponding upstream vari- 



ables are {xsu, cnu, fu) = (1.92, 



X 10" 



0.46). At large 



c, the numerical solution matches with asymptotic solution III8I I 
and ^ with A = 31.942 and B = 1.006 at large x. 



and the central thermal pressure decreases accordingly from 
10"^^ to 10"^^ dyne cm'^ 

We can also compar e vari able profiles with the case 
F of iTenorio-Tagle et al.l (Il986l ). The velocity profiles are 
very similar and we clearly see an expanding shock. As time 
evolves, the shock strength becomes weaker. In both nu- 
merical simulation and our model analysis, we observe a 
density peak on the downstream side of the shock and a sig- 
nificant temperature gradient on the upstream side, which 
cannot be accounted for by previous isothermal solutions. 
The upstream density and pressure profiles are also sim- 
ilar; however, the downstream density, pressure and tem- 
perature profiles (near ce ntre) are somewhat different. In 
iTenorio-Tagle et al.l ()l986l l. temperature, pressure and den- 
sity are initially very uniform behind the champagne shock 



but at the end of the calculations show large inward gra- 
dients. In Figure [T] our model also produces quasi-uniform 
temperature, pressure and density behind the shock at the 
begmnmg {t = 1.3 X 10^ yr), but we do not observe large in- 
ward gradients at the end. This differences are primarily due 
to the different physical assumptions adopted in the sim- 
ulatio n and our self-similar solutions. iTenorio-Tagle et al.l 
( 1986) treated the gas dynamics in protoglobular clusters 
and neglected the gas self-gravity as the gas mass is only 
about 0.1 per cent that of the stars. As shown above, our 
self-similar solutions neglect the gravity of the central mas- 
sive protostar but include the self-gravity effect. Another im- 
portant factor that introduces such difference is that in the 
simulation both forward champagne shock and reverse rar- 
efaction wave are taken into account. Our self-similar model 
can accommodate forward moving shock, so we only have 
the principle outgoing champagne shock. 

We calculate the total energy Etotai defined as the en- 
ergy of the gas under consideration in an infinite space. The 
total energy at time t can be expressed as 



-Etotai = Ek + Eg + El 



GMp 



+ |p liivr dr 



2G Jxi„ 



(3ri-2) 



a^x^{nx — v) 



dx 



(34) 



where Ek, Eg and Ei are the kinetic, gravitational and 
internal energies of the gas, respectively, and rin, rout, Xi-n, 
XoMt are the inner and outer boundaries of the gas under 
consideration, i is the degree of freedom of the gas parti- 
cle presumed to be 3. Note that one fixed rout at differ- 
ent times corresponds to different values of Xout. In this 
solution, Stotai = 5.5 x lO"** erg at t = 1.3 x 10^ yr, 
and iJtotai = 8.1 x 10'** erg aX t — 5.1 x 10® yr, so there 
is net energy input. Especially, we see the kinetic energy 
Ek = 6.4 X 10'*'' erg at t = 1.3 x 10^ yr, indicating a smaU 
fraction of the total energy, and Ek ~ 3.9 x 10*** erg at 
t = 5.1 X 10® yr, indicating a fairly large fraction of the 
total energy. The increase of kinetic energy shows clearly 
the development of a champagne flow. The gravitational en- 
ergy is of order 10''* erg, much less than the kinetic energy 
throughout this duration. This conflrms that the gas is not 
bounded and must have an outflow. We are also able to con- 
sider qualitatively the local energy exchange throughout the 
self-gravitating gas with relation (|16|l . In Figure [71 dp/dr is 
positive on the downstream side and negative on the up- 
stream side. For 7 = 1.1, the downstream and upstream 
sides locally loses and gains energy, respectively. In sum- 
mary, the proflles on the order of magnitudes, and the time 
evolution of our self-similar solution in modelling this case 
ar e grossly consistent with th e numerical simulation result 
of iTenorio-Tagle et al.l (|l986l ). which lends support to our 
polytropic self-similar "champagne flow" solution as a gross 
description of dynamics of H II regions. 

Recent studies further suggest that the inclusion of stel- 
lar winds is also important and even necessary sometimes 
in unders t andin g the large-scale dynamics of H II regions. 
IComeronI Ewh found that a shocked stellar wind in the 
central region produces important morphological differences 
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as compared to windless cases. Moreover, IComero 

3 (|l997l ) 

suggested that the spatial scale of an H II region undergo- 
ing "champagne flow" is systematically larger and the gas 
flow is generally faster as driven by a central stellar wind. 
[Arthur fc Roaxi (|2006t ) provided two-dimensional cylindri- 
cal radiative-hydrodynamic simulations of cometary H II re- 
gions using champagne flow models, by taking into account 
of strong stellar winds from the central ionizing star. In 
these simulations, the hydrodynamics and radiative trans- 
fer are coupled through an energy equation whose source 
term depends on the photoionization heating and radia- 
tive cooling rates; while with our polytropic approach, com- 
pli cated energetic process es are relegated to the choice of 
7. [Arthur fc Hoarg (|2006l ) studied the hydrodynamics of a 
compact H II region with a radius 0.13 pc and at a time 
~ 200 yr after the triggering of UV ionizing photons and 
powerful stellar winds; a stellar wind bubble around the 
centre with a radius up to 0.03 pc is formed. Inside such 
a stellar wind bubble, the mass density is about 3 orders of 
magnitude lower than that of the surrounding medium, and 
the density of the flow does not vary much with radius in the 
vicinity of the bubble boundary. Because the central wind 
bubble is effectively depleted of mass and the gravity force 
of the central massive star may be neglected, given a typical 
Bondi-Parker radius of ~ 10^ AU, we thus approximate such 
a stellar wind bubble as a central 'void' and model it using 
our polytropic self-similar void "champagne flow" model. 

In the scenario as outlined bv [Arthur fc Hoard (|2006l ). 
the central star has an effective temperature TcS = 3 x 10* 
K, a stellar wind mass-loss rate M = 10~® M0 yr~^ and a 
terminal wind speed = 2000 km s~^. The initial ambient 
medium has a number density of no = 6000 cm~"^ and a tem- 
perature of To = 300 K. In our self-similar model for "cham- 
pagne flows" , the radius of a void boundary is r* = k^^'^x*t". 
By taking n = 0.8, we have k^^^x* = 1.34 x 10^ cgs unit to 
obtain a r* = 0.03 pc central void at a time of t = 200 
yr. The n value depends on the energetic process, includ- 
ing plasma cooling and radiative heating, of the flow. We 
further estimate from relation (|33p for a downstream sound 
parameter kd ~ 2.5 x 10^^ cgs unit, and for a self-similar void 
boundary x* = 2.68. Another parameter that needs to be 
specified is the mass density on the expan ding void bound- 
ary, d enoted as a* here. The simulation of [Arthur fc Hoarg 
(|2006l l gives an electron number density on the void bound- 
ary as n* = 10* cm~^ at t = 200 yrs. With relation 
p* = a* /{AirGt^), we estimate a* — 5.5 x 10~^. We empha- 
size that the length and time scales in t his case are quite dif- 
ferent from those in the previous case of [Tenorio-Tagle et al.[ 
(|l986l l. As an important advantage, this suggests that self- 
similar models are suitable to give a unified description for 
cloud systems on quite different scales. 

As n < 1 in this case, we have one more degree of free- 
dom to specify the shock position. In principle, the shock 
position is determined by both the initial density (mass pa- 
rameter A) and the initial gas motion (velocity parameter 
B). We find that the lower limit of the downstream shock po- 
sition Xsd is 8.8, according to condition (|26|) : thus, the mini- 
mum shock position at t = 200 yrs is ~ 3.0 4 x 10^^ cm (i.e., 
about 0.1 pc). The numerical simulation of [Arthur fc Hoarg 
(|2006l ) studied the gas dynamics up to a radius of 0.13 pc. 
As an example of illustration, we assume x^d = 9 and show 
the resulting solution in Figure[8] This solution clearly shows 



that as time evolves, the void boundary expands, meanwhile 
the density and pressure in the vicinity of void boundary de- 
crease by several orders of magnitude, consistent with the 
simulation. However, we see a very high density near the 
"champagne flow" shock on the downstream side, and as 
time evolves, the density profile becomes more and more 
smooth. For our self-similar "champagne flow" model with 
central expanding voids, the velocity can rise up to several 
hundred km s~^ . In Figure[51 we also see clearly that the case 
is non-isothermal, and on the downstream side of the shock 
the temperature is the highest as expected. We note that at 
t = 800 yrs, the shock is at ~ 0.3 pc, beyond the scale of UC 
or compact H II regions. In reality, a champagne shock is so 
fast that even at a short timescale of ~ 800 yrs the shock is 
well in the surrounding diffuse interstellar medium (ISM). 

Compared with numerical simulations, the advantage 
of our semi-analytical self-similar approach is clear. We can 
generate self-similar shock solutions to model different H II 
regions by varying a few parameters. The self-similar pro- 
cesses and shock solutions of this paper describe the basic 
hydrodynamics of polytropic "champagne flows" and serve 
as test cases for bench marking numerical simulations. 

5.2 Asymptotic Free- Fall Solutions around a 
Central Protostar 

So far we have constructed "champagne flow" solutions with 
LP type asymptotic solutions on the downstream side as 
X ^ 0'^ , because such asymptotic solution satisfles bound- 
ary condition (|23p . Complementarily, free-fall asymptotic 
solutions (|20p and (|21|) at a:: 0^ represent gas infall 
and collapse during the protostar formation phase; and 
the surrounding gas and the infall momentum associated 
with the star formation process may be sustained for a 
while during the evolution after the onset of stellar nuclear 
burning and UV photoion ization of the surrounding gas. 
[Cochran fc Ostriked (|l977[ ) have investigated consequences 
of the birth of a massive star within a dense cloud with a 
free-fall density profile, and found that the radiation pres- 
sure from the star sweeps up grains from the infalling gas 
to form a dust shell which bounds the H II region. Here, 
we utilize such free-fall solutions as the downstream side 
and construct global solutions with shocks to model possi- 
ble dynamic evolutions of H II regions surrounding a nascent 
protostar in nuclear burning. We present such solutions in 
Figures [5] and [TD] where parameter m(0) for free-fall solu- 
tions is different. In dimensionless form, m(0) stands for a 
central mass point, and with dimensions in self-similar trans- 
formation dTj), M{0,t) cx f^""^m(0); therefore m(0) scales 
as the central mass accretion rate. For m(0) = 0.546 {Fig- 
ure[§]), the free-fall solution crosses the sonic critical curve 
smoothly ai x = 0.3237, and can also be connected to the 
upstream solutions via shocks at various locations. The free- 
fall solution crosses the line « = at a;stg ~ 0.74. This stag- 
nation radius expands with time in a self-similar manner; 
inside the stagnation radius the gas falls inwards, while out- 
side the stagnation radius the gas expands outwards. There- 
fore if Xsd < Xstg, the entire global solution corresponds to 
an inflow (solution 4 of Figure [9]). This situation describes 
an accretion shock during a protostar formation phase. If 
Sad > a^atg, the outer part of the downstream side is an out- 
flow. This describes the scenario that the shock sweeps up 
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Figure 8. Self-similar "champagne flow" solution with an initial central void radius of 10^'^ cm at time 200 yr (solid curve), 500 yr 
(dashed curve) and 800 yr (dotted curve). The four panels show mass density p, flow velocity u, thermal pressure p, and temperature 
T of the gas, respectively. The central void has a radius of x* = 2.68, corresponding to r* = 0.03, 0.06, 0.09 pc with increasing time t. 
The self-similar solution is obtained with parameters: n = 0.8, 7 = 1.2, a* = 5.5 X lO"'^ and downstream shock position Xsd = 9- The 
downstream sound parameter fed is 2.5 X 10^^ cgs unit, and the upstream sound parameter feu is 4.9 X 10^® cgs unit. The self-similar 
variables on the downstream side of the shock are {xg^, oj, v^) = (9, 1.07, 6.86), and the corresponding upstream variables are 
{xsu, Cfu, fu) = (20.4, 0.111, 8.88). At large x, the solution matches with asymptotic solution (1181 and (I19I I with A = 135 and B = 22 
at large x. The enclosed mass by radius of 10^^ cm is 1.38 X 10®, 1.23 X 10^ and 1.06 X 10® Mq as time evolves. The total energy of gas 
is 5.86 X lO'^^, 5.94 x lO^^^ and 5.98 x lO'^^ erg, respectively. 



the gas and turns the gas from infall to expansion on the 
dovirnstream side near the shoclt front. Similar to the situ- 
ation with downstream LP type solutions, there exists one 
specific Xsd, from which the upstream solution is a breeze 
with _B = 0. In this case, x^d = 1.7747 gives an upstream 
breeze (solution 2 of Figure (Qjl. Thus for Xsd < 1.7747, the 
upstream side corresponds to an asymptotic inflow far from 
the centre (solution 3 of Figure|5} and for Xsd > 1.7747, the 
upstream side corresponds to an asymptotic outflow (wind) 
far from the centre (solution 1 of Figure [9} . Another ex- 
ample shown in Figure \W\ has m(0) — 4.638. This free-fall 
solution does not cross the sonic critical curve smoothly and 
can be connected with upstream solutions via shocks. In an 
analogous manner, we show the possibility to obtain an out- 
flow (solution 1 of Figure llOp , breeze (solution 2 of Figure 
nop and inflow (solution 3 of Figure llOp for the upstream 
side. Solutions 1 and 2 of Figure [9] and solution 1 of Fig- 
ure [To] have asymptotic outflow or breeze on the upstream 
side and in the outer part of the downstream side, which 
is very similar to champagne flow solutions obtained with 
a downstream LP type, with different behaviours in central 
regions. With a free-fall asymptotic solution, the gravity of 
the central massive star is not neglected and the gas im- 
mediately surrounding the massive protostar still undergoes 
infall when the outer envelope starts to expand. Hence, solu- 
tions with free-fall centre are plausibly suitable to describe 



the early stage of "champagne flows". In general, with cen- 
tral free falls on the downstream side, we can also obtain 
asymptotic outflow, inflow, breeze and contraction for the 
upstream side, by varying the downstream shock position 
Xsd in a proper range. 

Central free-fall solutions describe the core collapse 
phase in star formation. For this purpose, IWang fc Loul 
( 2008) explored such solutions with free-fall inner core and 
an inflow or outflow in the outer envelope, for general poly- 
tropic cases in their Figure 2. The inner and outer portions 
are connected by magnetohydrodynamic (MHD) shocks. 
Such shocks are interpreted as accretion shocks, typically 
found in a star formation process, or around accreting black 
holes. Here we specifically emphasize that such shocks may 
also arise by the UV photoionization of ambient medium 
surrounding a nascent protostar. Under certain situations, 
the UV flux from the burning star might not be intense and 
rapid enough to turn the surrounding gas from infall to ex- 
pansion by ionization and heating. Meanwhile the ionization 
front (IF) creates a weak shock traveling outwards, and the 
upstream side may have an outward velocity. Another possi- 
bility is that the gravity of the central star is so large that the 
gas immediately surrounding the star keeps falling towards 
the protostar, but the outer part of the downstream side 
and the corresponding upstream side expand. In summary, 
with different initial conditions of gas and different physical 
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Figure 9. Reduced mass density a{x) (top) and reduced ra- 
dial flow velocity v{x) (bottom) for global solutions in cases 
with n = 0.9 (thus 7 = 1.1) whose downstream side is free- 
fall solution and the upstream side corresponds to either out- 
flow, inflow, breeze or contraction. In both panels, the dashed 
curve represents the sonic critical curve; in the bottom panel 
the dotted line is ii = 0. The downstream solution is connected 
with the upstream solutions with solid curves via shocks. The 
downstream solution is integrated from a sonic critical point 
(x, a, v) = (0.3237, 5.0050, -0.8455) towards x ^ 0+ with 
a central free-fall asymptotic solution of m(0) = 0.546, and out- 
wards to the downstream shock positions. At the inner most part 
the downstream solution corresponds to an inflow, and outer 
part of the downstream side is an outflow. The static point in 
this case is at Xgtatic ~ 0.74. In both panels, the upstream so- 
lutions from top to bottom correspond to x^^ = 2.5 (labeled 
1), 1.7747 (labeled 2), 1 (labeled 3) and 0.7 (labeled 4). Solu- 
tion 1 has inner inflow and outer outflow on the downstream 
side and upstream outflow with A = 4.368 and B = 1.76. The 
shock parameters are (x^di oj, v^) = (2.5, 0.8036, 1.3249), 
and {xsu, «u, Du) = (2.5003, 0.6453, 1.0981). Solution 2 has 
inner inflow and outer outflow on the downstream side and an 
upstream breeze with A = 1.8635 and B = at large x. The 
shock parameters are (x^^, a^, v^) = (1.7747, 1.0715, 0.8470) 
and {xau, Ou, Vu) = (1.7796, 0.5575, 0.1558). Solution 3 has 
an inner inflow and outer outflow on the downstream side and 
upstream inflow with A = 0.6930 and B = —1.6963 at large x. 
The shock parameters are (x^j, a^, v^) = (1, 1.7883, 0.2668) 
and (xsu, Ou, Vu) = (1.0118, 0.6348, -0.8942). Solution 4 has 
a downstream inflow and upstream inflow, with A = 0.4899 
and B = —2.1564 at large x. The shock parameters are 
(isdi Cfdi ''d) = (0.7, 2.4355, —0.0567) and (xsu, Qu, Vu) = 
(0.7054, 0.9837, -1.0785). 



conditions of a burning protostar, the radiative influence of 
the nascent protostar on the dynamic evolution of the sur- 
rounding gas may give rise to various self-similar solutions, 
including the classical champagne flow solutions, the ISECE 
solutions and the inner free-fall with outer inflow/outflows 
or contraction/breeze solutions. The solutions constructed 
in this paper (classical champagne flows) is suitable for sit- 
uations that the gas is initially static and the protostar ion- 




Figure 10. Reduced mass density a{x) (top) and reduced veloc- 
ity v{x) (bottom) for global solutions in cases of n = 0.9 (thus 
7 = 1.1) whose downstream side is a free fall and upstream side 
corresponds to either outflow, inflow, breeze or contraction. In 
both panels the dashed curve represents the sonic critical curve. 
The downstream solution is connected with upstream solutions 
of solid curves via shocks. The downstream solution is integrated 
from a sonic critical point (a; , a, v) = (1.7727, 1.0050, 0.5463) to- 
wards a: — » 0+ for a free-fall asymptotic solution of ni(0) = 4.638. 
Most of the downstream side is an outflow, while the inner most 
part is a free fall. In both panels, the upstream solutions from 
top to bottom correspond to x^^ = 1.6 (labeled 1), 1.4269 (la- 
beled 2) and 1 (labeled 3). The entire upstream solution la- 
beled 1 has an outflow with A = 2.8032 and B = 0.6058. 
Shock parameters are {x^^, a^, v^) = (1.6, 1.2326, 0.5110), and 
(xsu, «u, Vu) = (1.6002, 0.9575, 0.2443). Solution labeled 2 has 
an upstream contraction with A = 2.115 and B = at large x. 
Shock parameters are (x^d, "d, fd) = (1.4269, 1.5391, 0.4705), 
and (xsn, an, Vu) = (1.4290, 0.9057, -0.0988). Solution la- 
beled 3 has an inflow for the entire upstream portion with 
A = 1.0063 and B = —1.6727. Shock parameters are 
(Xjd, Odi '^d) = (1) 2.9346, 0.3423), and (xsu, Ou, Uu) = 
(1.0271, 0.7726, -1.2514). 



izes the entire gas immediately, and then the gas begins to 
expand in a "champagne phase" with an outgoing shock. 



6 CONCLUSIONS 

We present newly established self-similar polytropic shock 
solutions with and without central voids to model "cham- 
pagne flows" in H II regions featuring various asymptotic 
dynamic behaviours. As a substantial generalizat ion of the 
isothe rmal model of Tsai & Hsu (1995) and IShu et all 
we found similarities and differences in self-similar 
polytropic processes. Our general polytropic "champagne 
flow" model allows a much larger freedom to choose the 
polytropic index 7^1 for 2/3 < n < 2; for a conven- 
tional polytropic g subclass of examples, the power- 
law index / of the initial mass density profile p oc r~' is 
linked to 7 by i = 2/n — 2/(2 — 7). Together, our model 
is adaptable to a wide range of initial mass density profile 
with 1 < / < 3 for H II regions. For conventional polytropic 
cases of 1 < 7 < 4/3 (i.e., 2/3 < n < 1 and 2 < / < 3), we 
have more freedom for convergent initial conditions. In this 
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fashion, our conventional polytropic shock flow solutions are 
determined not only by the initial mass density profile (i.e., 
mass parameter A) , but also by the motion at the very early 
stage (i.e., velocity parameter B). The dimensionless shock 
positions or the dimensional shock speed and strength are 
determined by the initial conditions related to A and B and 
the central density at qq. Our self-similar shock flow solu- 
tions give a plausible description for the "champagne flow" 
phase for the dynamics of If II regions. We conclude that 
general polytropic "champagne flows" with the initial den- 
sity power-law index f < Z < 3 may evolve in a self-similar 
manner. 

We have established novel "champagne flow" shock so- 
lutions with an expanding void surrounding the centre to 
model a certain cloud core whose inner part has fallen into 
a nascent protostar. We observe that the evolution of central 
void boundary plays an important role in determining the 
asymptotic solution to approach and the general behaviour 
of solutions as well. With even one more free parameter, the 
"champagne flow" shock solutions with central voids can 
model the dynamics of H II regions more realistically, in- 
cluding the effect of central stellar wind bubbles. 

We have further explored possibilities of asymptotic in- 
flows or contractions far from the cloud centre. In addition, 
we also establish global shock solutions with the asymp- 
totic free-fall solution approaching the centre. In general, by 
varying dimensionless shock position, we connect the down- 
stream side, with either LP type solutions, EdS solutions 
or free-fall solutions, to upstream solutions which eventu- 
ally merge into asymptotic outflow, breeze, contraction and 
inflow. Within the theoretical framework of the self-similar 
polytropic fluid, global shock solutions with different be- 
haviours correspond to different forms of hydrodynamic evo- 
lution of H II regions after the nascence of a central mas- 
sive protostar. Apparently, even within the framework of 
self-similarity, dynamic evolution of polytropic H II regions 
depends on the initial and boundary conditions of molecu- 
lar clouds. Numerical simulations are needed to probe and 
connect various self-similar evolution phases. 
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APPENDIX A: AN INVARIANT FORM 

Transformation ((7)l leads to an initial density profile p oc 
r~^^" as the outer part of a SPS. To model an initial mass 
density profile other than a SIS (e.g., Galli et al. 1999) in 
the self-similarity framework, Shu et al. (2002) developed an 
invariant self-similar transformation for an isothermal gas by 
neglecting the self-gravity. We extend below this invariant 
transformation to a conventional polytropic gas. 

The independent similarity variable is a; = 
and the reduced radial velocity v(x) still reads as t; = 
M/(fc^/^r-^). But the mass density p is generalized to 

p(r,t) = ^i?(a;) , (Al) 

where K{x) is a new reduced density dependent on x only, 
/ is a scaling index and D is a constant parameter to be 
determined. With the conventional polytropic relation p = 
Kp^ for constant k and 7, the thermal pressure p is 

p(r,t) = K^W . {A2) 
We now determine the relation between k and k from the 
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du ^ du 
dt dr 



(A3) 



radial momentum equation without self-gravity 

Idv 
p dr 

The left-hand side (LHS) of equation (|A3|) relates to the 
sound scaling parameter k and the right-hand side (RHS) 
relates to the specific entropy coefficient We substitute 
self-similar transformation (lAip and (IA2|I into both sides of 
the radial momentum equation (IA3|) to obtain 



^+u^= k^^h"-^ [(n - l)v + (v- nx)v'] , (A4) 
dt or 



Idp 
p dr 



nD-'-^-f{k-^^H-'') 



l/2,-n\{7-l)i + l -(7-l)i 



Lr-<-^ + R'R-<'^ ] , (A5) 



where the superscript '/' over v{x) and R{x) indicates the 
first derivative in x. By expressions HA4|) and (|A5|) and in 
order to remove the explicit t dependence in the self-similar 
form of equation (|A3|I . we require 

- n[(7 - 1)Z + 1] = n - 2 . (A6) 

Together with n -|- 7 = 2 for a conventional polytropic gas, 
we obtain I = 2/n for n ^ 1. For n = 1, there is no con- 
straint on I. For a conventional polytropic gas, the possible 
initial mass density profile described by this invariant self- 
similar transformation is not arbitrary, but is preset by the 
scaling index n or by the polytropic index 7. Here D pro- 
vides a measure for the magnitude of density profile and 
K is the constant in the polytropic equation of state; both 
parameters are determined physically. For simplicity of the 
formulation, we may choose the sound parameter k so that 

n = . (A7) 

This relation is verified by a dimensional analysis. As a re- 
sult, these parameters disappear in the following ODEs, 



= —\a - 2)iv - nx) + (71 - l)x] - 'ylx^-'R'' 

X 



(A8) 



where H is an arbitrary constant but should be chosen for a 
positive mass density everywhere. The condition I — 2/n has 
already been taken into account in this solution. This form 
of asymptotic solution at large x differs from equation (34) 
of Shu et al. (2002), because for non- isothermal cases the 
leading terms in nonlinear ODEs (|A8|I and (|A9|I are differ- 
ent; however, the second term on the RHS of the expression 
of V, reduces to the v expression in equation (34) of Shu et 
al. (2002) for n = 1 and 7 = 1. The form (|1I0)) is simi- 
lar to the form of (|18|) and (|19p in the overall scaling. The 
integration constant H in the invariant form is an analogy 
of the 'velocity parameter' B, while the 'mass parameter' 
A is already absorbed into the parameter D in the invari- 
ant form. To make this asymptotic solution applicable, we 
should require H = in cases of n < 1 and thus 7 > 1. 
The boundary conditions at x ^ -|-oo correspond to the 
initial conditions of the fluid. Therefore, the initial density 
scales as r~' . More precisely, this new self-similar transfor- 
mation requires the initial density profile to scale as r~^^" , 
identical with the scaling we obtain in the main text using 
self-similar transformation (0. For x ^ 0^, the asymptotic 
solution reads 



' 2x/3 



R ~ Ix' 



0' 



(All) 



where / > is an arbitrary parameter for a positive mass 
density. This v expression of asymptotic solution (|A11[I is 
independent of I, different from the v expression (35) in 
IShu et all (|2002D . but it is similar to asymptotic solution 
H24|) in the main text; the R expression of our asymptotic so - 
lution is similar to the R expression (35) of lShu et al.l ||2002D . 
and the parameter I here is proportional to the central re- 
duced mass density ao. In the isothermal case with n = 1 
and 7 = 1, relation (|A6|I is satisfied automatically and pa- 
rameter I is thus arbitrary. Here lies a fundamental difference 
between the conventional polytropic case and the isothermal 



[(nx — n)^ — 72;^ 'i?'' ^]v' 

= 2'yx^~' R'^~^ {v — x) + {n — l)v(nx — v) 



(A9) 



For n 7^ 1 and I = 2/n, coupled equations (IA8|I and HA9|) 
describe non-isothermal self-similar flows of a conventional 
polytropic gas. For n = 1 and 7 = 1, the value of I is arbi- 
trary; for I = 2, these two coupled nonlinear ODEs reduce to 
the isothermal formulation of Shu et al. (2002) with two de- 
coupled ODEs for I — 2. For n — 1 and 7 = 1 with Z 7^ 2, we 
need to use the two decoupled ODEs (32) and (33) in Sec- 
tion 4 of Shu et al. (2002). In contrast to the isothermal case, 
the nonlinear ODEs here remain coupled for 7 7^ 1. In sum- 
mary, with our self-similar transformation, hydrodynamics 
for a spherical gas with ignorable gravity can be expressed 
by a set of two coupled nonlinear ODEs (|A8|I and (IA9p . 

In cases with 7 7^ 1 and for x — » +00, the asymptotic 
boundary conditions are v ^ and R 1 for n < 1. The 
asymptotic solution at large x yields 



+ ■ 



R^1 + 3H{ - 

n 



1 X 



-l/r. 



+ 00 , 



(AlO) 



